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ABSTRACT 

In  a  cylindrical  region  -we  consider^electromagnetic  fields  independent  of 
the  axial  coordinate:  controlling  the  time  evolution  of  such  fields  by  means 
of  boundary  currents,  likewise  independent  of  the  axial  direction,  is 
equivalent  to  controlling,  simultaneously,  two  wave  equations)  one  with 
boundary  control  of  Dirichlet  type,  the  other  of  Neumann  type.  In^his  paper 
QjmT provides^  preliminary  study  of  control  problems  of  this  type  and  indicate 
what  is  necessary  for  extensions  of  our  work. 
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SIGNIFICANCE  AND  EXPLANATION 

This  paper  concerns  the  controllability  of  the  Maxwell  electromagnetic 
equations  in  a  cylindrical  spatial  region  by  means  of  controlling  currents 
caused  to  flow  on  the  boundary  of  the  region.  Here  controllability  refers  to 
the  ability  to  transfer  from  electric  and  magnetic  fields,  given  at  the 
initial  instant,  to  corresponding  fields  prescribed  at  a  later  instant. 

Studies  of  this  type  are  significant  in  relation  to  wave  guides,  EM-pulse 
devices,  radar  non-relective  (stealth)  aircraft,  controlled  thermonuclear 
fusion  and  many  other  important  applications. 
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THE  DIRICHLET-NEUMANN  BOUNDARY  CONTROL  PROBLEM  ASSOCIATED 
WITH  MAXWELL'S  EQUATIONS  IN  A  CYLINDRICAL  REGION 


D.  L.  Russell 

1.  BACKGROUND. 

In  this  paper  we  consider  a  region  (1  Q  R3,  not  necessarily  bounded,  having  piecewise 
smooth  boundary  T  and  almost  everywhere  uniquely  defined  unit  exterior  normal  vector 
v  •  v(x,y,z),  (x,y,s)  e  T.  It  is  assumed  that  the  region  a  is  occupied  by  a  medium 

having  constant  electrical  permitivity  e  and  constant  magnetic  permeability  u .  We  have 
then,  in  0,  the  paired  electric  and  magnetic  fields 

E  -  E(x,y, z,t)  , 

H  »  H(x,y,z,t)  , 

having  finite  energy 

*(t)  -Vi/J/  (elil2  +  jilSl2 )dv  ,  <1.1) 

a 

where  I  ■  denotes  the  usual  Euclidean  norm  in  R3.  As  is  well  known  ([4],  [9]),  £  and 

H  satisfy,  in  Q,  Maxwell's  equations 


♦ 

Be 

curl  H  “ 

e  Jt  ' 

(1.2) 

-  ♦ 

3h 

curl  E  “ 

(1.3) 

div  £  - 

P  , 

(1.4) 

div  H  - 

0  , 

(1.5) 

where  p  *  p(x,y,z,t)  is  the  electrical  charge  density  in  a  -  which  is  zero  throughout 
this  paper.  (That  aquation  (1.5)  might  eventually  have  to  be  modified  to  account  for 
magnetic  monopoles  will  trouble  us  not  at  all  here  I ) 
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Control  problms  associated  with  Maxwell's  equations  have  been  of  interest  prinarily 
in  connection  with  nuclear  fusion  applications  -  in  which  case  p  is  not  identically  equal 
to  sero  and  the  Maxwell  equations  are  coupled  with  the  dynamical  equations  governing  the 
plasma  evolution.  In  this  connection  we  cite  the  work  of  P.  K.  C.  Kang  [29],  [30],  [31]. 
The  point  of  view  which  we  take  here  is  that  we  cannot  hope  to  treat  these  more  complicated 
problems  until  we  have  a  firmer  grasp  on  the  control  theory  of  Maxwell's  system  in  its  own 
right.  In  this  direction  some  work  on  controllability  with  control  Influence  distributed 
throughout  Q  has  been  carried  out  by  G.  Chen  [2],  [3].  We  are  primarily  concerned  here 
with  the  possibility  of  influencing  the  evolution  of  the  fields  £  and  fi  by  means  of  an 
externally  determined  current  5( x,y , z, t )  flowing  tangentially  in  T  so  that 

$(x,y,s,t)*3(x,y,z)  -  0  ,  (1.6) 

for  (x,y,z)  e  r  where  v(x,y,z)  is  defined.  We  will  assume  that  the  normal  component  of 
i  vanishes  outside  fl  and  that  no  charge  is  permitted  to  accumulate  on  T .  Then  we  have 
the  boundary  conditions  (see  e.g.  [4],  [28]) 


c£(x,y,z,t)*v(x,y,z)  »  0  (1.7) 

U$T(x,y,z,t)  “  v(x,y,z)  *  J(x,y,z,t)  (1.8) 

for  (x,y,z)  e  T  such  that  v(x,y,z)  is  well-defined.  Here,  and  subsequently,  the 
subscript  T  refers  to  the  component  of  the  vector  in  question  which  is  tangential  to  V. 
Similarly,  the  subscript  v  will  denote  the  normal  component  (thus  (1.7)  is  the  same  as 
"  0).  writing 


£*■£  +  f  »£  on  r, 

V  T  I 

fi  -  fiv  +  fit 

1  »  $T  on  r  , 

we  see  that  (1.8)  becomes  pST  *  V  x  so  that  is  a  vector  tangential  to  r  and 

perpendicular  to  3  m  3  • 

The  state  space  in  which  we  study  solutions  of  the  above  system  will  be  denoted  by 
d(0)i  it  is  a  closed  subspace  of  the  space  HE(0)  of  square  integrable  six- 
dimensional  fields  (£(x,y,z,t),  S(x,y,z,t))  with  the  inner  product  and  norm 
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<1.9) 


<<l1,it1)»<l2,iJ2)>  =  ///  <e£,*l2  +  u^itjJdv 
l(S,fi)l2  -  <(l,fi);(ljfi)>  . 


Clearly  Hg(ft)  is  a  real  Hilbert  space  with  this  inner  product.  Where  a  complex  space  is 
required,  we  employ  conjugation  as  usual.  The  state  space  d<0)  is  the  closed  span  in 
Hg(fl)  of  those  continuously  differentiable  fields  (E(x,y,z,t),  H(x,y,z,t))  for  which 


div  E 


div  H 


SE 


3e 


3E_ 


TT  *  W  +  TT 


3H 


3H 


3H 


TT  +  w  *  TT 


o  , 

o  . 


If  £0.ff0  and 
the  first  corresponding  to 


are  two  smooth  solution  pairs  for  { 1.2)  —  ( 1.5),  (1.7),  (1.8), 
J  =  0  on  T,  we  see  easily  that 


ft  <(^o'"o,;<*i'"i,>  “ 

si,  si-  si,  3l 

"  III  (etEo*  W  +  It-  ’  EJ  +  ft-  +  3t“  *  "J)dv 


(using  (1.2),  (1.3)) 


“IS!  <V  curl  H  i  “  curl  E0  *  +  curl  1^*  1^  -  Ip*  curl  l^dv 

“  (using  div  (E  x  H)  ■  curl  I  •  I  -  I  •  curl  I) 


-  ///  tdiv(!0  x  f^)  +  divtl,,  x  l0)]dv 


-  //  <SQ  x  H,  +  l1  x  H0)  •  vds  -  (using  (1.7)) 


-  //  <S0T  -  *lT  -  e0t  «  «1v  ♦  Stt  *  l0T  ♦  l1t  *  !0v,  .  vd. 
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-  {/  <S0T  *  *  S1t  *  V  *  Vd* 


(using  (1.8)  and  noting  that  J  S  0  for  2..  S. ) 


//  (I  •  S)ds 

r 


(1.10) 


If  we  go  through  the  same  computation  with  ^g'^g’^i'^i  both  replaced  by  the 
satisfying  (1.2)-(1.5),  (1.7),  (1.8)  we  find  that 


S  ■  -  JJ  (I  *  f)  •  vda  "  -  JJ  2  •  Jds  .  (1.11) 

«  r  r 

For  J  S  0  generalised  solutions  of  (1.2)-(1.S),  (1.7),  (1.8)  can  be  discussed  in  the 
general  context  of  partial  differential  equations  and  strongly  continuous  semigroups.  The 


generator 


A(l,2)  »  (j-  curl  2,  -  ~  curl  I) 


with  domain  consisting  of  2.2  in  the  Sobolev  space  n  1,1  having 


sero  divergence  and  satisfying  (cf.  (1.7),  (1.8)) 


*.lr  -  ”•  Mr  ■  “  • 


is  antisymawtric  and  generates  a  group  of  isometries  in  H^d(fl).  (See  f 32 1 ,  [33],  [34] 
for  related  work.)  Sufficient  conditions  on  2  so  that  solutions  of  the  inhomogeneous 
system  (1.2)-(1.5),  (1.7),  (1.8)  lie  in  d((J)  and  are  strongly  continuous  there  may  be 
obtained  much  as  in  [18],  [19]  but  it  is  not  easy  to  specify  necessary  and  sufficient 
conditions.  Indeed,  this  is  already  difficult  for  the  much  sinq>ler,  but  related,  wave 


equation 


a2w  _  a2w  t  a2w  i  a2w 
3t2  3x2  3y2  3s2 


with  boundary  forcing  terms.  He  will  make  some  comments  related  to  this  in  Section  6. 
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2.  CONTROL  PROBLEMS  IN  A  CYLINDRICAL  REGION 

The  main  point  in  thia  paper  ia  to  study  the  question  of  controllability  of  the 
electromagnetic  field  by  means  of  the  boundary  current  3  m  3^.  By  controllability 

we  mean  the  possibility  of  transferring  an  initial  field  l(x,y,z,0),  &(x,y,z,0)  e 
HB,d'°)'  9iven  at  time  t  •  0,  to  a  prescribed  terminal  field  E(x,y,z,T), 
fitx,y,z,T)  e  Hg  ^(Q) ,  specified  at  t  »  T  >  0,  by  means  of  a  suitable  control  current 
3(x,y,z,t)  defined  for  (x,y,z)  e  V,  t  e  (0,1).  Because  the  homogeneous  Maxwell 
equations  correspond  to  a  group  of  isometries  in  Hg  d<0),  it  is  enough  to  consider  the 
special  case  wherein 

t(x,y,z,0)  =  0  ,  (2.1) 

S(x,y,z,0)  i  0  .  (2.2) 

ror  a  given  space,  J,  of  admissible  control  currents  3(x,y,z,t)  »  ^(x,y,z,t)  defined 
on  T  x  [0,T]  we  define  the  reachable  set  JUT, J)  to  be  the  subspace  of  Hg  ^(Q ) 
consisting  of  states  reachable  from  the  zero  Initial  state  using  controls  3  e  J. 
Following  earlier  definitions  <(8],  (26)),  our  system  is  approximately  controllable  in 
time  T  if  R(T,J)  is  dense  in  Hg  d(Q)  and  exactly  controllable  in  time  T  if 
R(T,J)  «  hb  d<0)  (or  some  precisely  designated  subspace  of  Hg  d(B)). 

At  this  writing  we  are  not  able  to  discuss  the  general  three  dimensional  problem 
wherein  the  vector  fields  £  and  3  are  unrestricted,  except  as  stipulated  heretofore, 
and  Q  has  a  general  geometry.  Ne  hope  in  later  work  to  consider  at  least  some  three 
dimensional  cases  which  arise  for  special  domains  Q.  But  for  now  we  must  content 
ourselves  with  the  case  in  which  0  is  a  cylinder: 

0  ”  R  *  (-*•,«)  ■  {(x,y,z)  |  (x,y)  e  R  c  R2,  z  real) 
where  R  is  an  open  connected  region  in  R  with  piecewise  smooth  boundary  B.  Thus 

30  -  3R  x  (-»,*•)  -  B  x  (-»,»)  . 

Even  here  we  can  give  results  only  for  special  two  dimensional  regions  R. 

The  two  dimensional  problem  in  the  cylinder  0  «  R  x  (-»,»)  occurs  when  we  confine 
attention  to  fields 

I  ”  E(x,y,t),  H  ”  H(x,y,t) 
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which  do  not  depend  on  tha  coordinate  t  corresponding  to  the  axial,  or  longitudinal, 
direction  of  tha  cylinder.  (Note  that  this  is  not  at  all  the  sane  thing  as  requiring 
that  *  ,  H  ,  the  field  components  in  the  z  direction,  should  be  zero. )  We 
correspondingly  consider  only  control  currents 

5  »  5(x,y,t) 

which  do  not  depend  upon  z. 

Of  course  the  energy  S  in  Q  is  infinite  under  the  above  circumstances  if  t,  9 
are  not  identically  zero.  We  redefine  I  to  be  the  energy  per  unit  length  of  cylinder: 

I(t)  "^JJ  (elf(x,y,t)l2  ♦  gl9(x,y,t)l‘ )dxdy  .  (2.3) 


The  space  d(fl)  is  now  replaced  by  Because 


3E  (x,y ,t)  3H  (x,y,t) 

H  °'  ^  5  0 


we  have 


3E  3e  +  3H  3H 

div  ®  “  J5T  +  5y  '  dlw  H  “  3^  +  By*  * 


The  curl  expressions  simplify  to 

curl  E  3y  '  5 3x  3y  '  ' 

3Hz  3^  3H 

curl  H  -  (sj-,  -  YT ’  3T  '  W>  ' 
so  that  the  equations  (1.2),  (1.3)  become 


<  1 )  c 


(ii)  C 


3E 

3H 

3H 

x 

z 

(iv) 

X 

U  - — 

atT  " 

3y 

3t 

3E 

3H 

3H 

— y . 

at 

z 

”  3x 

(v) 

“3^ 

3E 

3H 

3H 

3H 

z 

■Ml  V 

_JC  . 

X 

(vi) 

z 

p  r— 

3t 

3x 

3y 

3t 

!!!• .  -  !!z  ♦  !!a 

3t  "  3x  3y 


It  is  clear  from  (2.5),  (i)-(vi),  that  if  f(x,y,0),  H(x,y,0)  are  given,  then  the 
subsequent  evolution  of  Ez(x,y,t),  Hz(x,y,t)  determine  all  of  the  other  components.  As 
for  these  components  themselves,  differentiating  (2.5)  (iii)  and  (2.5)  (vi)  with  respect 


to  t  and  then  substituting  (2.5)  (iv),  (v)  and  (2.5)  (i),  (ii)  into  the  respectively 
resulting  expressions,  we  obtain  the  familiar  wave  equations 


32e 

z 

32  Ez 

32E 
►  _ £ 

3t2 

3x2 

3y2 

32H 

z 

32H 

„  _ £.  H 

32H 
,  _ * 

3t2 

3x2 

3y2 

valid  for  (x,y)  e  R,  t  e  [<),«•),  provided  Ez,  Hz  have  enough  derivatives,  or  provided 
the  equations  are  interpreted  in  the  distributional  sense.  Assuming  the  initial  states 
i(x,y,0),  H(x,y,0)  are  divergence-free,  we  compute  (cf.  (2.4)) 


and  similarly 


»  Ob  Ob 

E  9t  +  3^)  “  (using  (2.5)  (i),  (ii)) 

32H  32H 

”  3^^  “  0 


3_ 

3t 


3H  3H 

(jr +  ^ 


and  we  conclude  that  the  fields  remain  divergence-free  for  all  time. 

Suppose,  then,  that  divergence-free  initial  states  S(x,y,0),  H(x,y,0)  are  given. 

3E 

Then  Ez(x,y,0),  Hz(x,y,0)  are  known  and  (2.5)  (iii),  (vi)  determine  (x,y,0)  and 


3H 
_ 1 

3t 


(x,y,0).  If  (2.6),  (2.7)  are  then  solved  with  these  initial  conditions,  and 


appropriate  boundary  conditions,  the  complete  solution  of  Maxwell's  equations  (2.5) 
(i)-(vi),  can  he  obtained  by  integrating  (2.5)  (i),  (ii),  (iv),  (v).  Thus  it  is  enouc^i  to 
work  with  (2.6),  (2.7),  and  it  should  be  noted  that  the  divergence  condition  does  not  have 
any  bearing  on  Ez,  Hz»  it  can  be  ignored  henceforth. 


nn  J.n. :» 


Figure  1.  The  Region  ® 


It  ia  important  to  recast  the  boundary  conditions  (1.7),  (1.8)  so  that  they  provide 
boundary  conditions  for  (2.6),  (2.7).  We  ask  the  reader  to  consult  Figure  1,  where  the 
region  R  with  boundary  3R  *  B  is  shown.  At  a  point  (x,y)  e  B  we  let  $  “  v(x,y) 
denote  the  unit  exterior  normal  to  B  and  we  let  o  “  o(x,y)  denote  the  positively 
oriented  unit  tangent  vector  to  B  there.  With  the  unit  vector  in  the  positive  z 

direction,  v,  $,  (  form  a  positively  oriented  orthogonal  triple  of  unit  vectors.  Given 
an  arbitrary  vector  w  we  can  decompose  it  as 

*  ■  (vvV-V>  ' 

,♦,2  _  2  .  2  .  2 

Iwl  "  wv  ♦  "g  *  • 

The  tangential  part  of  ft.  which  we  have  designated  as  ftT>  may  now  be  represented  as 

ftT  -  +  H^o  (2.8) 

and  the  current  1  "  $T  may  likewise  be  represented  as 

3x  “  Js*  +  ' 
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(2.9) 


vkJ-v*Jt=ux  <J^J  +  J(jS)  -  -J^S  +  Jgt 
Combining  (1.8),  (2.8),  (2.9)  we  see  that  on  B 

Hz(x,y,t)  -  J0(x,y,t)  , 

H0(x,y,t)  -  -J^tx.y.t)  . 

♦  + 

Represent  v,  a  as 


Then  compute 


v  *  v  t  *  v  n  , 
x  y 

3  =  o?  +  ori  =  -v$  +  vri. 

X  y  y  X 


3e 

z 

3v 


3e 


3e 


- — ~  v  +  — —  v  »  (using  (1.3)  , 
3x  x  3y  y 


3H 


3H 


y  .  x 

u  a  +  y  *—  <j 

3t  y  3t  X 


!! *s 

3t 


=  (using  (2.11)) 


3J 

z 

3t 


(2.10) 

(2.11) 


(2.12) 

(2.13) 


(2.13) 


(2.14) 


The  equations  (2.10),  (2.14)  provide  the  needed  boundary  conditions  for  (2.6),  (2.7) 
respectively.  For  Hz  we  have  the  Dirichlet-type  boundary  condition  (2.10)  while  for 
Ez  we  have  the  Neumann-type  boundary  condition  (2.14).  If  we  let 

U(x,y,t)  =  (x,y,t)  , 

5  =  6  =  (I  5  M)  l  , 

t  <y  z 

and  differentiate  (2.10),  we  have  the  more  symmetric  form 


3h  3e 

(x,y,t)  =  U0<x,y,t),  «  -lMx,y,t),  (x,y)  e  B  .  (2.15) 


We  complete  this  section  by  discussing  the  question  of  expression  of  the  energy  per 
unit  cylinder  length,  (2.3),  solely  in  terms  of  Hz  and  Ez> 

we  consider  the  equations  (2.6),  (2.7)  with  homogeneous  boundary  conditions 

3h  3e 

tt  <x'y't)  “  °'  rr  <x-y<t)  ‘  (x'y)  e  B  • 
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We  um  the  symbol  4  for  the  Laplacian: 


Initially  we  take  Hz,  Ez  to  lie  in  the  Sobolev  space  H2(R).  This  space  Bust  be 
decomposed  in  order  to  attach  a  meaning  to  i” 1 . 

The  boundary  condition  for  Hz  may  be  rewritten  as 

Hz(*,y,t)  "  h(x,y),  (x,y)  e  B  , 

where,  by  the  trace  theorem,  h  e  H3/,2(B).  Then  we  can  write 

Hz(x,y,t)  »  H^tx.y.t)  +  H^(x,y ) 
where  H^(x,y)  is  the  solution  of 

AHz<x,y)  »  0,  H^tx.y)  «  h(x,y),  (x,y)  e  B 

and 

H^(x,y,t)  »  0,  (x,y)  e  B  . 

—  1  * 

The  inverse  Laplacian  4  is  well  defined  on  the  functions  H^.  For  Bz  we  may  write 

Ez(x,y,t)  »  Ez(x,y,t)  +  E^lt) 

where  E^,  as  indicated,  is  constant  with  respect  to  (x,y)  e  ft  and 

/  E  (x,y,t)ds  *  0  . 

B  * 

-1  "* 

It  is  well  known  that  A  is  well  defined  on  the  functions  E  . 

z 

We  proceed  first  on  the  assumption  that 

H,(x,y,t)  *  H  ( x,y ,t 9 ,  E_(x,y,t>  =  E  (x,y,t)  . 

4  Z  z  z 

We  form  new  solutions  of  (2.6),  (2.7)  by  setting 


3G 

z 

w 


-E 


3Fr 

3t 


G 

z 


eA 


.  3E 
-1  _ z 

3t  ' 


UA 


.-1  ^s 
3t 


We  then  determine  Gy,  Gy,  Fx,  Fy,  using  the  equations  (2.5)  with  £  replacing  B,  t 
replacing  £.  so  that  P  and  £  satisfy  Maxwell's  equations: 


-curl  F 


It  will  then  be  found  that 


3f  * 

c  *  curl  G  . 


E  »  curl  F,  H  "  curl  £  . 


Following  this,  (2.3)  can  be  written  as 


*(t)  »  V2  /  J  (elcurl  $l2  +  ulcurl  &l2)dxdy 

s 

3F  2  3F  2  3F  3f  2 

*(5T)  *  <)/  '  5Tl  1 


3G  2  3G  2  3G  3G  2 

*  '«**)  *  (af)  ]d*ay 


Then  from  (2.16)  we  have 


3f  2  3F  2  3G  2 

*(t)  *V2//  {€[(3^)  +(^)  Mm  z~)  1 

3G  2  3G_  2  3F.  2 

+  m[(— )  +  Cay-3  +  (e  jr)  H  **r 


3F  2  3F  2  , 

h/M<sr>  ‘(if)  ♦<v1i 

R 


3G  2  3G  2 


+  +  ^3T}  +  ‘V  ]dXdY 


Now  consider  the  quadratic  form  (for  Ez  =  E^) 
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Similarly 


it  ,  it  32G  ,  J2G 

(IT*  “  it)  •  <■  rr*  -*'1  -r) 


“  (lince  G  satisfies  the  wave  equation  lie  — 5-  *  4G 

3t^  * 

and  the  boundary  conditions  Gz(x,y,t)  »  0,  (x,y)  e  B)  - 

,  ,  3G  2  3G,  2 

—  (-40  .g  )  — j  Un-J  +  (jf)  1  * 

»e  ue  7 


3H  .  3h  ,  3r  2  3F  2 

(IT*  IT)  •  t  Wit)  ♦  W  1 


from  which  it  follows  that 


-  3E_  ,  3B_  3H  f  3H 

,2rr  *  *■>  .  t  * 


*(t)  “Vi/J  {(i»e)  [(yp>  “4  yj^)  ♦  (jt-*  yjr)]  +  etEB'  +  M(Hi)  i****1* 


I  I  finite  states  -  a  fact  which  will  be  very  careful  later. 

It  is  necessary  to  aodify  this  expression  for  general  Ez,  Hz.  We  begin  with 

B  tx,y,t)  “  l  (t)  . 

9E* 

The  only  possible  solutions  of  the  wave  equation  (2.6)  satisfying  j^~|B  “  0  and  having 
this  fora  are 

K,(x.y.t)  -  e0  +  e,t 

where  eQ  and  e^  are  constants.  (Such  solutions  are  consistent  with  a  constant  boundary 
current  J  for  which  Jg  =  0 . )  The  corresponding  tx.  By,  Hz  are  rero  but 

3E  3h  3h 
t  y  x 

*i  •  It-  5x  "  ?F  * 


It  is  not  possible  to  express  this  quantity  in  terns  of  Ez  itself  or  Hz.  It  is  better 

3E 

z 

to  leave  It  in  the  form  £  xr— •  Solution*  of  Maxwell's  equations  with  Ez  having  this 

9E 

z 

fora  have  energy  expressible  as  a  quadratic  form  in  E  and  r— • 
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■ 


N«xt  we  consider  Hz  •  as  described  earlier.  Such  a  solution  is  consistent  with  a 
boundary  current  for  which  JT  ■  0,  constant  with  respect  to  tiM  but  possibly  varying 
with  (x,y)  e  B.  We  may  take  Hx,  Hy,  Et  all  sero.  However, 


3Ex  3HX  3E 

e  W  ■  ST'  6  it 


3x 


so  we  may  not  assume  that  Ex  and  Ey  are  equal  to  zero.  The  energy  associated  with 

solutions  of  this  type  is  expressible  in  terms  of 

3H  2  3H_  2 

II  [(f^-)  +  (fT^)  ]«Sxdy 


if  integration  with  respect  to  t  is  permitted.  In  the  sequel  we  will  not  explicitly 
consider  the  timewise  linear  electric  fields  satisfying  the  above  equations. 

We  see  then  that  a  norm  involving  only  Ez  and  Hz  and  compatible  with  the  energy 
(2.3)  may  be  expressed  as 


dE  9e  3h  3h 

<<VV'2  “  //  -A  1  |^)  +  -A  1  ♦  c<*z)2  +  W(Hr)2 

R 


VV'  +  ’1&2  +  *  •IIS5)' W 


(2.17) 


where  p0,p 1 ,a0,a,  are  positive  numbers.  It  will  be  seen  that  this  is  a  weaker  norm  than 
the  one  associated  with  a  pair  of  wave  equations,  viz.: 

9E  2  3h  2 

,(VHz)|2  *  //  Mtft^  +  (ft5)  ]  *  ,VKr|2  +  'VHtl2}dxdy  .  (2.18) 

3E  3h 

We  will  denote  the  Hilbert  space  of  states  Ez,  Hx,  j—i  lying  in  h’ir),  h’ir), 

2  2  * 

L  (R),  l‘(r),  respectively,  by  H.  This  space  will  be  very  convenient  for  uee  in  the 


remainder  of  this  paper.  In  some  cases  we  will  add  boundary  conditions  to  the 
specification  of  H,  the  space  with  norm  I  I,  without  changing  the  symbol,  to  correspond 
to  an  agreed  specification  of  the  states  in  H  by  similar  boundary  conditions. 


\ 

► 


3.  SOW  CONTROL  CONFIGURATIONS 

*•  describe  here  two  possible  realization*  of  th*  control  problem  which  we  have  posed 
and  indicate  why  we  have  chosen  the  mathematically  more  interesting  (i.e.,  more  difficult) 
one  to  work  with  in  this  paper. 

Let  us  assume  that  r  *  )C  •  l  «  (-«•,••) 
is  covered  by  one  or  more  layers  of  conducting 
bars,  arranged  in  rows  as  shown  in  Figure  3.1. 

In  the  case  of  a  single  layer  of  conducting 
bars  shown  in  Figure  2(b),  the  bars  are  arranged 
so  that  they  make  an  angle  9,  0  <  |8|  <  i, 

with  the  vector  o  (cf.  Figure  1),  while  in 
the  double  layer  case  (Figure  2(a))  they  are 
arranged  so  that  the  bars  in  the  second  layer 
make  an  angle  0,  0  <  |0|  <  0  *  8,  with 

the  vector  I.  The  current  in  any  row  of  bars 
parallel  to  the  z-axis  is  independent  of  z; 
i.e.,  constant  for  all  bars  in  that  row.  As 
we  consider  successively  smaller  bars  we 
obtain,  as  an  idealization,  the  boundary 
current  vector 

S(x,y,t)  »  J(x,y,t) (cos  00  +  sin  0^)  (3.1) 

in  the  single  layer  case,  J(x,y,t)  denoting  the  current  strength  with  the  sign  determined 
so  that  J  positive  yields  a  positive  current  component  in  the  o  direction.  The 
corresponding  formula  in  the  double  layer  case  is 

3(x,y,t)  »  J.j(x,y,t)(cos  00  +  sin  0c) 

+  Jj(x,y,t)(cos  00  +  sin  0c)  .  (3.2) 

The  current  components  are,  in  the  single  layer  case 

Ja<x.y.t)  ■  J(x,y,t)cos  0  , 

J^lx.y.t)  “  J(x,y,t)sin  0  , 
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Figure  2(b).  Single  Layer  Control 


j*. .. 


1 


I 

I 


-ww 


and  In  tha  double  layer  case, 

J0<x,y,t) 

Jjtx.y.t) 


coa  8  coa  t 
aln  8  sin  V 


J^x.y.t) 
J2(x,y,t)  . 


(3.3) 


The  determinant  of  the  matrix  In  (3.3)  Is  sin  (^>  -  8)  4  0  if  If  *  8  In  the  range 
0  <  |8|  <  §,  0  <  Itl  <§.  Thus  in  the  double  layer  case  Jg  and  are  independent 

if  J>|  and  J2  are  independent  while  in  the  single  layer  case  Jg  and  Jz  are  fixed 
non-zero  multiples  of  each  other. 

The  double  layer  case  is  easily  disposed  of  in  the  light  of  earlier  work  on  boundary 
control  of  the  wave  equation.  Referring  back  to  (2.10),  (2.11)  we  now  have,  for 
(x,y)  e  B  -  9R,  t  e  (0,»), 


3H 

j—  (x,y,t)  -  Oc(x,y,t)  “  cos  8  ut(x,y,t)  cos  +  u2<x,y,t)  , 


ff-  (x,y , t )  -  -U^  (x,y,t)  -  -sin  8  u1(x,y,t)  +  cos  <>  u2(x,y,t)  , 

3J,  3J2 

ut(x,y,t)  -  (x,y,t) ,  Uj(x,y,t)  «  j~-  (x,y,t)  . 


Since  Vg  and  Uz  are  independent  if  uj  and  u2  are,  the  control  problem  splits  into 
two  uncoupled  wave-equation  problems,  one  for  Ez  and  one  for  H^.  These  have  been 
discussed  thoroughly  in  [2],  (3],  [15],  [16],  [22],  [23],  [25]  with  affirmative 
controllability  results  for  various  control  configurations  and  will  not  concern  us  further 
here. 

In  the  remainder  of  this  paper  we  study  the  single  layer  case.  If  we  let 

u(x.y.t)  -  (x.y.t) 

we  now  have  the  wave  equations  (2.6),  (2.7)  for  E£,  Hr  and  the  boundary  conditions 


(3.4) 


3v 


(3.6) 


.1  m _ i  - 


3j 


The  control  problem  for  >c  and  H(  ara  now  coupled  bacauaa  tha  single  control 
function,  u(x,y,t),  appears  in  tha  boundary  conditions  for  both  «8  and  Kci  we  have  to 
contol  both  system  siaultanaously  using  tha  earn  control  function. 

If  wa  rely  on  experience  in  a  single  space  dimnsion,  which  has  proved  generally  quite 
helpful  in  the  control  theory  of  a  single  wave  aquation,  we  are  led  to  believe  that  system 
like  (3.6),  (2.7),  (3.5),  (3.6)  my,  in  fact,  be  controllable.  Replacing  u(x,y,t)  by 
Ug(t),  Uf(t)  and  taking  0  <  x  <  1,  the  one  dimensional  equations  are,  using  variables 


v,  w. 


0 


(3.7) 


|^T  (0,t)  “  «uQ(t),  (1,t)  *  au1(t)  , 


(3.8) 


o  — 2 - 2  ■  0 

3t2  3x‘ 


(3.9) 


( 0,t)  -  -6u0(t),  ||  (1,t>  -  $u,(t) 


(3.10) 


(note  that  -  corresponds  to  tha  exterior  norml  derivative  at  0).  Letting 


~  Jv 
3x 


(3.11) 


we  find  that 


and 


"*  3w 
w  -  — 
3t 


,  2~  ,2- 

.  3  v  3  v  B 

0  — 2 - 5  “  0  ' 

3t  3x* 


„2~  .2~ 

3  w  3w  . 

0  — 2 - 2  “  ° 

3t  3xZ 


(3.12) 


(3.13) 


(3.14) 


Differentiating  (3.11)  with  respect  to  t  and  using  (3.8)  we  have 
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.7  * 

■  : 


-  (0,t>  -  -  (O.t)  -  -  g'lt)  , 

P  ,x2  P  9*  P  0 


2i_S  (1<t)  «  2  |Z  (1t)  .  •  „<(t|  . 

P  ,.2  P  9x  pi 


while  differentiation  of  (3.12)  along  with  (3.10)  yields 


(0*t)  ■  §f  (O.t)  -  -Bu'(t)  . 


-£d,t>  -Bu;(t, 


Coebining  (3.13)  with  (3.14),  (3. IS),  (3.16),  (3.17),  (3.16),  we  see  that 
Bv  +  2  w<  Bv-^w  both  satisfy  the  wave  equation  and 


(Bv  *  2  S)( O.t)  -  0,  ~  (8v  +  S  5)(i,t)  -  u.(1) 

3x  p  3x  p  pi 


(3.16) 


(3.17) 


(3.16) 


f;  (Bv  -  2  5)(0,t)  -  uj(t),  g  (Bv  -  2  7X1.0  -  0  . 


Thus  the  control  problems  for  Bv  +  —  w  and  Bv  -  —  w  are  both  of  Kausa nn  type  and  are 
uncoupled.  Affirmative  controllability  results  are  then  available  froai  [201,  [21],  [24]. 

If  we  replace  u^lt)  (or  u1(t)  by  0  in  the  above,  then  Bv  -  2  w  (or  Bv  +  2  w) 
will  becoae  completely  uncontrollable  and  our  original  system  must  therefore  be 
uncontrollable.  This  result  at  first  seems  to  predict  failure  for  the  enterprise  which  we 
now  undertake  for  the  two  dimensional  case. 


4.  APPROX IHAT1  BOUNDARY  CONTROLLABILITY 

By  •  simple  change  of  acala  in  tha  t  variable,  and  renaming  of  the  independent 
variables,  we  say  assume  that  the  ay  a  tea  of  intereat  ia 


3t2 

ifw 

at2 


afv  +  afv 

ax2  ay2  ' 


a2«  a2w 
ax2  3y2 


t  >  0  , 
(x,y>  e  R  , 


(4.1) 


(4.2) 


with  boundary  conditiona 


(x,y,t )  -  au(x,y,t) 
|~  (x,y,t)  -  8u(x,y,t) 


t  >  0  , 

<x,y)  e  B  -  3Q 


(4.3) 


(4.4) 


We  will  not,  in  general,  assume  that  u(x,y,t)  can  be  aelected  at  will  for  all  valuea  of 
(x,y,t)  ahown.  More  on  thia  later. 

Becauae  the  system  ia  time  reversible,  it  ia  sufficient  to  analyse  controllability  in 
terms  of  control  from  the  sero  initial  state 


v(x,y,0)  -  (x,y,0)  -  0 


w(x,y,0)  -  (x,y,0)  »  0  , 


V  (x,y)  e  R  > 


(4.5) 


(4.6) 


to  a  final  state 


3v 


v(x,y,T)  ■=  Vg(x,y ) ,  (x,y,T)  «  v^x.y) 


aw 


w(x,y,T)  -  wQ(x,y) ,  ^  (x,y,T)  -  »,(x,y) 


(x,y>  e  R  . 


(4.7) 


(4.8) 
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Me  have  noted  In  Section  2  that  the  I  I-  finite  statea  are  dense  in  the  I  I-  finite 


states.  In  the  present  context  this  means  that  we  can  work  with  the  Hilbert  space  of 

3v  3w 

states  v,  w,  —  with  the  inner  product 


(( 


v. 


3v 

It' 


3  v 

It' 


■//  [ 


3v  3v 

Jt  It 


+ 


3w  3w  ^ 

It  It  + 


3v  3 v 

K  IS 


(4.9) 


a  space  which  we  will  refer  to  as  H.  The  norm  is  I  I  (cf.  (2.18))  with  lit  ■  1.  As  we 
have  indicated,  this  is  a  dense  subspace  of  H,  the  Hilbert  apace  obtained  by  use  of  the 
norm  I  I  (cf.  (2.17)). 

The  final  states  (4.7),  (4.8)  are  not  quite  arbitrary  in  H  if  the  control  u  is 
restricted  so  that  its  support  is  contained  in  a  proper  relatively  closed  subset  C  B. 
Since  the  condition 

(x,y,t)  ”  a  u(x,y,t),  (x,y)  e  B 

applies,  we  may  as  well  adjoin  the  additional  condition 

v„(x,y)  -  0,  (x,y )  e  B  -  B,  5  Bq  .  (4.10) 

The  trace  theorem  ([11,  [19])  assures  us  that  this  describes  a  closed  subspace  of  H, 
which  we  will  call  H^.  The  only  restriction  on  H^  is  (4.10) i  v0  is  permitted  to  have 
arbitrary  values  in  H,/'2(b1)  and  wg,  w.,  are  unrestricted  in  H^B),  H°(B)  “  L^(B), 

respectively. 

Let  U  be  a  qiven  space  of  admissible  control  functions,  about  which  we  will  shortly 

have  more  to  say.  For  each  control  u  6  U  we  assume  the  existence  of  a  unique  solution 

Vu,  wu  of  (4. 1 )— (4.6)  for  t  »  0,  (x,y)  8  B.  Very  general  sufficient  conditions  for  this 

to  be  the  case  are  given  in  [19].  We  define  the  reachable  set  at  time  T,  R(U,T),  to  be 

3v  3w 

u  _ u 

the  set  of  all  final  states  vu(x,y,T),  (x,y,T),  wu(x,y,T),  (x,y,T)  which  may  be 

realized  in  this  way.  The  set  R(U,T)  is  a  subspace  of  Hf  if  U  is  a  linear  space, 
which  we  will  assume,  and  our  system  is  approximately  controllable  in  time  T  if  R(U,T) 

is  dense  in  H^  (then  R(U,T)  is  also  dense  in  H  because  II  is  a  weaker  norm  than 

•  * 

I  I  and  is  dense  in  H).  Evidently  R(U,T)  is  dense  in  H^  just  in  case,  given  an 


arbitrary  atata  (*o'*1,*0'“l>  in  H1' 

3  vu  3wu  ^  ^  ^ 

i(C*u<**y»T>'  (x.y.T)/  **u(*.y,T),  ~  (x,y,T));(vQ,v1,w0,w1))  -  0  , 

u  e  u)  —  >  <VVV-1>  “  °  ■  (4.11) 

Let  v(x,y,t),  w(x,y,t)  ba  tha  uniqua  solution  of  (4.1),  (4,2)  satisfying  the  terminal 
conditions  at  time  T» 

*-  *»•  3v  *0  •»  3w 

v(x,y,T)  -  vQ,  y^  (x,y,T)  -  v1#  w(x,y,T)  -  w0,  y^  (x,y,T)  -  v»1  ,  (4.12) 

and  the  hoaegenaous  boundary  conditions 


<x,y,t)  -  0  , 
|~  <x.y,t)  -  0  , 


(x,y)  e  B,  t  >  0  . 


(4.13) 


(4.14) 


Cosiputing  tha  quantity 


«  ((v**y.t>.  w1  (x'y't)'  wu(x«y't}'  »r  1 

(v(x,y,t),  y“  (x,y,t) ,  w(x,y, t ) ,  (x,y,t>))  , 

using  familiar  duality  theorems  involving  tha  Laplacian  and  integrating  from  0  to  T 
(sea  (22],  [23],  [26]  for  details  in  the  case  of  a  single  wave  equation)  we  see  that 


3  v  3w 

((vu(x,y,T>,  y^  (x,y,T)  ,  Wu(x,y,T),  y—  (x,y,T))>  <v0  #V,  »w0  >**,  >) 


T  «**  3v  a~  3v 

/  /  [57  (x,y,t)  — *  (x,y,t)  +  57  (x,y,t)  y~  (x,y,t) 


0  B 


+  I7  (x,y,t)  y-2  (x,y, t)  +  (x,y,t)  (x,y,t)]dsdt 


3w-  3v  3w 


at 


aw 


av 


at 


(4.15) 


Then  using  the  boundary  conditions  (4.3),  (4.4),  (4.13),  (4.14)  we  see  that  the  above 


-20- 


raducaa  to 


Following  the  development  In  [6) ,  it  nay  be  lean  that  our  ay  a  tarn  is  exactly 
•*  2 

controllable  in  B, ,  uaing  the  control  apace  0  -  L  (B^  x  [0,Tl),  juat  in  case 


+  >  kkv  v 

8  8t  L Z(B1xtO,Tl)  H 


(4.20) 


for  some  K  >  0.  In  general  this  is  a  very  difficult  result  to  obtain  but  we  are  able  to 

2 

obtain  exact  controllability,  by  other  means,  for  the  case  where  B  is  a  disc  in  R 
and  B1  ■  B  is  its  boundary,  a  circle.  This  result  is  developed  in  -iction  6  where  it 
will  be  seen  that  it  is  heavily  dependent  on  certain  properties  of  the  Bessel  functions. 
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The  work  here  can  be  carried  out  for  a  rectangle  with  arbitrary  dimensions,  but  all 


essential  ideas  are  contained  in  the  notationally  simpler  case 

*  =  { (x,y ) | 0  <x<*,  0<y<*} 

to  which  attention  is  restricted  henceforth.  We  will  assume  that  B.) ,  the  portion  of  the 
boundary  on  which  control  is  exercised,  is  one  side  of  R,  without  loss  of  generality  it 
is  the  set 

B,  =  { (*  ,y)  1 0  <  y  <  w)  .  (5.1) 

We  consider  then  v,  w  satisfying  (4.1),  (4.2)  in  R  *  [0,T)  for  some  T  >  0,  and  al3o 
satisfying  boundary  conditions 

( x,y, t )  =  0,  (x,y,t)  =  0,  <x,y)  e  B  =  3R,  t  e  [0,T]  ,  (5.2) 

“  1Z  +  8  If  I''*'" 

=  “  lx  +  8  If  “°»  0  <  y  <  t  e  [0,T]  .  (5.3) 


We  may  assume  without  loss  of  generality,  since  the  wave  equation  is  time  reversible 
with  either  Dirichlet  or  Neumann  boundary  conditions,  that  v  and  w  are  extended  to 
satisfy  (4.1),  (4.2)  on  -*»  <  t  <  "  and  that  the  boundary  conditions  (5.2)  hold  for 
(x,y)  e  B,  t  e  (-*>,<»).  We  may  not  assume  that  the  boundary  condition  (5.3)  is  applicable 
beyond  [0,T],  however,  if  controls  are  restricted  to  have  support  in  B1  x  [0,T).  Let 
5  >  0  and  let  s(t)  be  an  arbitrary  function  in  C  (-*»,«•)  with  support  in  (-5,5). 


Define 


v(x,y,t)  »  /  s(t  -  T)v(x,y,T)dT  , 

-OS 

A  » 

w(x,y,t)  «  /  s(t  -  T)w(x,y,i)dT  . 


(5.4) 


(5.5) 


Then  v,  w  are  solutions  of  the  wave  equations  (4.1),  (4.2)  satisfying  boundary  conditions 


dv  d  w 

—  (x.y.t)  -  0,  —  (x.y.t)  •  0,  fx,y)  e  B  •  3B,  -«•  <  t  <  •  ,  (5.6) 


a  (x.y.t)  +  $  (x.y.t)  -0,  0  <  y  <  x,  t  (  U,  T  -  !)  .  (5.7) 


Moreover,  it  can  be  shown  that  v,  w  are  of  class  C  for  (x,y)  e  R»  -*•  <  t  <  “.  If  we 
can  show  v  =  0,  w  =  0  for  any  such  choice  of  s,  then  v  s  0,  w  =  0. 

Let  us  define,  for  (x,y)  e  R,  -*•  <  t  <  », 


From  (5.7)  we  have 


♦(x.y.t)  «  a  (x.y.t)  +  0  (x,y,t)  . 


♦(x.y.t)  *0,  0  <  y  <  x,  t  e  [6,  T  -  8] 


Since  a  and  6  are  constants  we  have 


=  +  (x,y)eR,  _<t< 
3t  3x^  3 y* 


Let  us  note  that,  since  v  satisfies  the  wave  equation  in  RUB, 


a  i-j  (x.y.t)  *  8  (x,y,t ) 

3t 


2a  2^  2* 

(x.y.t)  +  (x.y.t)]  +  8  (x.y.t) 

3x  3y 


Setting  x  *  x  in  (5.11)  and  differentiating  the  identities  in  (5.6)  with  respect  to  t, 
we  see  that  the  left  hand  side  vanishes.  Then,  comparing  (5.11)  with  (5.8) 


f£  (x.y.t) 


-a  (x.y.t)  5  a(y) ,  0  <  y  <  x,  5  <  t  <  T  -  «  ,  (5.12) 

3y 


the  last  identity  being  valid  as  a  consequence  of  the  first  condition  in  (5.6). 

The  two  conditions,  (5.8)  and  (5.12),  satisfied  by  ♦  at  the  boundary  x  »  x  enable 


us  to  use  Holmgren's  uniqueness  theorem  (see  [5]  or  (13],  e.g.)  in  much  the  same  way  as  it 
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r 


! 


was  used  in  the  proof  of  the  approximate  controllability  of  the  wave  equation  in 
to  aee  that  if 


[22], 


[23] 


T  >  2  +  26  (5.13) 

then  ^  muet  be  independent  of  t  for  1+3<t<T-1-4,  i.e. 

♦(x,y,t)  -  ♦(x.y),  (x,y )  e  R,  1  ♦  J  <  t  <  I  *  1  -  (5.14) 

Because  v  and  w  satisfy  the  wave  equation  in  it  with  the  homogeneous  boundary 
conditions  (5.6),  and  are  of  class  C  in  R  U  B,  we  have  C*-  convergent  expansions 


v(x,y,t)  » 


w(x,y,t) 


a  ,  b  w.  ,  v. 

v0(x,y)  +  l  l  <v  e  3  +  v  e  3  )sin  kx  sinjy 

k-1  j-1  3  “3 

.  “  ltd  t  -  -i“k1t 

”  wn  +  I  l  (wvte  +  3  )c°8  k*  COB  iy  ' 

k=*1  j»1  3  ■’ 


(5.15) 


(5.16) 


where 


wk3  =  /7T7 , 


vq (x,y )  is  a  C  function  in  RUB  such  that  (cf.  (4.10)) 

Vq (x,y)  =  0,  (x,y)  e  B  -  {(*,y)|0  <  y  <  «> 

and  wQ  is  a  constant.  Then,  from  (5.8), 


(5.17) 


(5.18) 


♦(x,y,t)  -  a 


3vfl(x,y) 


3x 


l  cos  kx[  l  (akvk:)sin  jy  +  iB^^^cos  jy)e 


•  -  -iw.,t 

+  l  (akv^sin  jy  -  iBw^^^cos  jy)e  3  ]  , 


(5.19 


still  C  -  convergent  for  (x,y)  e  R  U  B,  ■*•  <  t  <  •.  Noting  (5.14),  we  see  that  the 
left  hand  side  takes  the  form 

4  A 

3vQ(x,y)  3vQ(x,y) 

*(x,y,t)  -  a  — - *(x,y)  -  o  - ^ -  S  ♦(x.y)  . 

1+4<t«T-1-4.  (5.20) 
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H«  now  strengthen  (5.t3)  to 


(5.19)  gives 


8v0(x,y) 

♦<x,y,t)  »  ♦(x.y)  «  a  - - ,  (x.y)  e  R  , 


(5.25) 


1+S<t<T-1-«. 

Noting  (5.15)  and  (5.16)  and  the  fact  that  v(0,y,t)  5  0,  we  conclude  from  (5.23)  that 


v(x,y,t)  5  vQ(x,y)  , 


w(x.y,t)  2  w 


0  ' 


1+«<t«T-1-4. 


(5.26) 


Since  v(x,y,t)  2  vQ(x,y)  is  a  solution  of  the  wave  equation  with  (cf.  (5.18)) 
vQ(x,y)  -  0,  (x.y)  e  B  -  {<»,y)|0  <  y  <  *} 

it  must  In  fact  be  a  solution  of  Laplace's  equation  there.  Then  we  cosqpute 


2* 

3v„ 


82v„ 


3v  2  3v  2  w  w  » _ 

I  ((^-  (x,y))  +  (3—  (X.y))  +  v0(x,y)(— (x.y)  ♦  — y  (x,y))]dxdy 

R  y  3x  3y 


/  div(vQ(x,y)grad  vQ(x,y))dxdy 


*  -  *  *  8v0 
/  vg  (x.y)  grad  v0(x,y)*V(x,y)dS  *  /  v0<w,y)  (*,y)dy  . 


(5.27) 


Combining  (5.9)  and  (5.25)  with  the  fact  that  v^  satisfies  Laplace's  equation  we  conclude 
from(5.27)  that 


3v  2  3v.  2 

/  [(f^-  +  (3^-  (x.y))  jdxdy  =  0 


and  this,  together  with  (5.18),  implies 


vQ(x,y)  2  0  . 


(5.28) 


Combining  (5.26)  and  (5.28)  we  conclude  that 


v(x,y,t)  2  0 
w(x,y,t)  2  w 


(x.y)  e  R,  -  -  <  t  <  - 


(5.29) 


the  result  for 


<  t  <  » 


being  an  immediate  consequence  of  the  result  for 


6.  SOME  EXACT  CONTROLLABILITY  RESULTS  IN  THE  CASE  OF  A  CIRCULAR  CYLINDER 


We  consider  now  the  case  fl  ■  R  *  ( -m,m)  with 

R  -{(x,y)|x2  +  y2  <  1}  i 
B  -  3R  »  {<x,y)|x2  *  y2  -  1)  . 

with  introduction  of  the  usual  polar  coordinates  r ,  6,  the  equations  (4.1),  (4.2)  now 
become 

•!:  .  !?i  ♦  2  |s  ♦  i- 1!=  (6.,, 

at2  3r2  r  3r  r2  ae2 


a2w  _  a2w  |  i  »v  t  i  a2w 
at2  3r2  r  3r  r2  ae2 


(6.2) 


and  the  boundary  conditions  (4.3),  (4.4)  are  transformed  to 


( 1,0, t)  -  au(8,t)  , 


3w 


( 1,6,t)  -  0u(8,t)  . 


(6.3) 

(6.4) 


Writing 


v(r,0,t)  -  l  vk(r,t)eikB, 


V-k  -  Vk  ' 


(6.5) 


#(r,8,t)  -  l  wk<r,t)eike. 


W-k  "  "k  ' 


(6.6) 


u(8,t)  “  J  u^ttje 


ik0 


(6.7) 


we  arrive  at  an  infinite  collection  of  control  problems  in  the  single  space  dimension,  r: 


*\  3\  1  3vk  k2  „  ,  k  ,  . 

TT  m  TT  4  x  Tz - 2  vk  ■  °'  -<*<-' 

at  3r  r 

3S  3\  1  3wk  k2  „  ^  ,  k  ,  . 

.2  .2  r  Sr  2  k 

at  3r  r 


(6.8) 


(6.9) 
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* 


M.t>  «  au^lt ), 

3wk 

—  d.t)  -  Su^t), 


—  <  k  <  •  , 


-  •  <  k  <  •  . 


(6.10) 


(6.11) 


We  will  first  treat  the  equation  (4.1)  with 
have  seen,  reduces  to  the  set  of  problems  (6.8), 


the  boundary  condition  (4.3)  which, 
(6. 10),  -•  <  k  <  •.  With 


as  we 


z(r ,6,t) 


zfc(r,t)e 


ik8 


3vk(r>t)  ik8 
3t  * 


3v 

3t 


(r ,6, t ) 


we  have  the  equivalent  first  order  systems 


3_ 

3t 


(r.t) 

(r,t) 


) 


IlfvJr.th  _  .  fv.  (r.t), 
0Jl*Nr,t)J  **|k|  ^zNr.t)' 


(6.12) 


where  L|  k  |  is  the  differential  operator  on  the  right  hand  side  of  (6.8).  The  boundary 
conditions  (6.10)  become 

zk<1,t)  -  cnik(t),  -»  <  k  <  •  .  (6.13) 

The  eigenvalues  of  the  operator  I.|k|  with  the  corresponding  homogeneous  boundary 
condition 

zk(1,t)  “  0  (6.14) 

are 

°'  ±i4i|k|,t'  *  “  1'2'3"--  ' 

where  “|k|  (  is  the  t-th  positive  zero  of  the  Bessel  function  J|k|(r)  of  order  |k|. 
The  corresponding  vector  eigenfunctions  are 


where 


♦|k|o(r»,  ,  ♦|kU(r)  , 

0  tl“lk,,t  *|kU(t,J' 


—  <  k  <  «• 
t  -  1,2,3,...  , 


♦|k,,o(r) 


Sk|,0r 


lk| 


—  <  k  <  «• 


(6.15) 
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4|k|t(r)  "  A|k|tJ|k| *W|k| ,lr^' 


-“  <  k  <  •  , 
t  -  1,2,3,.. 


The  normalization  coefficients  A|k|,0'  *  |  k  |  (  are  chosen  so  that 

I  r|*|k|0(r)|2dr  -  /  r  |  ♦ ,  k ,  #  j(r  > !  2dr  -  t-  1,2,3,. 


..  .  (6.16) 


Thus 


A  -  /in.  ± J 

|k|,0  /  w  ' 

while,  as  may  be  seen  from  (5),  e.g. 


-  «  <  k  <  *  , 


(6.17) 


Ikl.t 


^  ,ik|,*(**|k|,t> 


(6.18) 


The  state  space  in  which  we  wish  to  work,  for  the  present  at  least,  is  (cf.  (2.18)) 

H  -  {^)|v  e  h\*),  z  e  l2(»)J 

with  the  inner  product 

iO).  (i2))  *  /  ♦  z ) dxdy 

1  2  R 

and  associated  norm.  Since  the  f|k|  j  satisfy  the  homogeneous  boundary  condition  (6.14) 
one  easily  sees  that 

.  ik0  ,  _ 

|(  ,klo°e  )L  -  -  /  ♦|K|,0eik9A1*|k|.0ei,'8)d^ 

H  R 


.  r  .  _ik8  3*|k 1 ,0  _-ik6 .  |k|  +  1  ,2* 

31,  ,k|'°  5  q  Ik,dfl 


-  2 1 k  |  ( | k |  +  1),  —  <  k  <  -  , 


(6.19) 


while 
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1 


tor  T  >  0  we  arrive  at  the  equations 


ip  2w  d  ^ 

2 1  k  |  ( |k|  +  1)v„  .<T)  -  a  /  /  u(8,t)  -i-k1,0  (1)  e"ik9d8dt 

K#u  0  0  *r 


3*|fc1,o  T 


2*a  — (i)  /  uk(t)dt 


(6.23) 


21 


Ik| 


T  2v 

,v*  ,<T)  -  a  /  /  u(8,t  )e 

»*  K'*  0  0 


i“|k|,t(T-t>  8*|k|,i  ...  -ik8 


3r 


(1)  e  d8dt 


2va 


J*lkl  1  ,T  iUlkl  l(T_t) 
— (1)  /  e  ,k|'1  v 


(t)dt  , 


(6.24) 


2X|kl,*V;.i<T)  “«/  f  U(6't> 


T  r21 _ -i“|k|,t<‘r't>  ‘♦ikl.t  -ik8 


3r 


(1)  e  d8dt 


3r 


(6.25) 


Thue  the  Dirichlet  boundary  control  problem  for  (6.8),  (6.10)  1b  reduced  to  a  moment 
problem  (6.23),  (6.24),  (6.25)  for  which  ^(t)  must  be  a  solution.  We  proceed  in  much 
the  same  way  with  the  Neumann  boundary  control  problem  for  (6.9),  (6.11).  We  let 


C(r,8,t) 


Ck(r,t)e 


ik8 


3wk(r,t) 

3t 


ik8 

e 


3w 

3t 


(r,8,t) 


and  obtain,  in  place  of  (6.12), 


3 

It 


fwk 


(r,t) 

(r,t) 


) 


I 

0 


<r,t) 

(r,t) 


)  “  " 


lk| 


(r,t) 

(r,t) 


)  • 


(6.26) 


The  boundary  conditions  are  now 


3wk 

(l.t)  -  8^ (t) ,  —  <  k  <  -  . 


The  eigenvalues  of  M| ^ |  with  the  corresponding  homogeneous  boundary  condition 
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are,  for  k  •  0 


<1.t>  -  0 


0,  ±lv. 


t  ■  1,2,3,...  , 


where  Vg  j  is  the  t-th  zero  of  the  differentiated  Bessel  function,  jg(r),  of 


0,  and,  for  k  *  0, 


*iV|k|,t*  *  “  1'2'3'"'  ' 


where  t  is  the  t-th  zero  of  J^(r ) .  In  the  case  k  «  0  the  eigenvalue  0  has  double 
multiplicity.  The  special  solutions  taking  the  place  of  (6.22)  in  this  case  are 


w(-,*,t) 


C(*,*,t) 


<t  -  Tl*. 


where  is  such  that  (cf.  (6.161) 


, 1  .2  „  1  1 

I  r*oo  *  “  27'  A-e*  *oa  '  ~ 


In  all  of  the  other  cases  the  vector  eigenfunctions  take  the  form 


+Ik|.t(r) 


±lv|k|,t*|k|,t 


(r))«  —  <  k  <  »,  t  -  1,2,3,... 


<'|k|.t<r)  '  B|k|.tJ|k|(V|k|,tr)'  t-  1,2,3,.'..  , 


the  normalization  coefficients 


selected  so  that 


,  ,2,  1/2  ,  . 
**  (U|k|,t  ‘  k  ’  J|k|<V|kL* 


J  r|*|k|<t(r)|2dr  -  ^  . 


The  corresponding  special  solutions  of  the  homogeneous  equation  are 


2|1 


T  2*  iv,  ,  ,<T-t> 

|k|,twk,t<T)  “  6  l  l  ut9't)iv|k|,ie  '  *M,t 


(1)e‘lk9d8dt 


_  T  iv,  ,  ,(T-t) 

2l,W*|k|.«*tM.*‘1)  I  e  '  Vt,dt 


T  2» 


2“|k|,t<.l(T)--6/  l  u(9't)iV|k|  .£« 


'iu|k|,£<T“t)  me'ik9d0dt 
*|kl,t(1)e  d9dt 


_  T  -iv.  .  (T-t) 

~2,,Biv|k|,**|k|.£(1)  {  e  ’  Vt)dt 


We  find  also,  taking  (~)  in  the  second  form  given  in  (6.27),  that 


T  2*  _  _  T 

C00(T>  »  8  /  /  u<8,t)*00d8dt  -  2«$*00  /  u0(t)dt  . 

0  0  0 


Since  this  must  be  true  for  all  T  and  gr  woo^^  *  *  we  ^ave  a*so 


«00(T)  -  2W8<C00  /  (T  -  t)u0(t)dt 


Since  U(k,  £  »  (v,k,  £)  .  (6.31),  (6.32)  become 


»,.  .  ,  .  _  T  iv..,  ,<T-t) 


‘ikiyikl.t^lkl.t’  l  e  Ik1'1  Uk<t,dt 


n,.|  ,  T  -iv.  .  |  (T-t) 

> 1  v  (T)  •  -B,  .  J.  ,  .(v.,.  .)  /  e  '  u.(t)dt 

»8i  k,£  lk|  ,t  |k| ,£  |k) ,£  q  k 


Taking  account  of  the  fact  that 


*♦ 


Mil 


9J 


9r 


M>  -  <■>,,.  ,A 


Mil 


|k| ,£  |k| ,£  8r  '“|k|,t 


(6.24)  and  (6.2S)  yield 


(6.31) 


(6.32) 


(6.33) 


(6.34) 


(6.35) 


(6.36 
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U  (t)dt 
k 


(6.37) 


L  v+  (T)  .  A  _ lk|.d  (M  )  /T  . 

*o  k,l'T’  |k|,t  3r  |k|,i’  JQ 


T  i“|k|.t(T"t) 


"(kl  £  -  9jlvl  t  T  _i“lvl  ,(T_t) 

V  (T>  -  A  — (u  )  /  e  |k,,t  u  (t)dt 

*o  k,l  |k|,£  3r  |k|  ,t  '  k 


On  the  other  hand 


ao  (6.23)  gives 


- 1" (  1 )  -  A  Ikl 

3r  '  '  "|k | ,0  * 


Ikl  ♦  1 
wo 


vk.o<T)  -  A|k|,o  l  Vt,dt 


Using  the  formula  (6.18)  and  (6.28)  for  A j ^ |  t  and  B|k|  t  we  have 


W|fcl.*  w+  . ,T)  „  Vlkl.t  fT  iV|k| ,l<T‘t> 

wBi  k,i  ft  (M|k|  t  -  k2)V2  0 


ufc(t)dt 


-l*!/.*:  w“  (T)  -  - 

-Bi  k,Z(  ft 


.  -  *  >  '  o 


;T  e"iV|k| 


u  (t)dt 
k 


1*M  v+ 

wo  k 


.  T  i».  .  (T-t) 

(T)  “  — -  /  e  u  (t)dt 

.*  JZ  n  k 


ISli*  v'  <T)  -  —  / 

wo  k,£  y-  ' 


e  uk(t)dt 


The  equations  (6.39)  become,  in  view  of  (6.17), 


Ikl  (JkJ  +  1) 


\.o(T>  ’  ^  l  Vt)dt 


(6.40) 


This  is  valid,  but  meaningless,  for  k  -  0.  It  is  easy  to  see  that  in  the  case  k  »  0  we 
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should  use 


/w 


00 


(T)  -  /  u  ( t )  dt 
0  k 


(6.45) 


The  equations  (6.33)  and  (6.34)  are  left  as  they  appear.  We  note  that  all  of  the 
coefficients 


1X1,1 


_  _1  /2|k| 

r  2  1/2'  r-'  r- 


>  k  *  0,  2»$ 


(6.45) 


are  bounded  away  from  zero,  uniformly  with  respect  to  k. 

It  is  also  possible  to  show,  using  the  work  [10],  [11]  of  K.  D.  Graham,  that  the 
numbers 

°'  U|k|,1'  “lkl.1'  V I  k  |  ,2  '  “lkl.2 . V|k|,j'  “|k|,j"" 

are  separated  by  a  gap  at  least  equal  to  »/2  again  uniformly  with  respect  to  k. 

Applying  the  result  [14]  of  A.  E.  Ingham  along  with  the  work  of  Duff in  and  Schaeffer  [7], 
much  as  in  [12],  [2],  [3],  we  conclude  the  existence  of  functions  uk(t)  in  L2[0,T],  for 
any  fixed  T  >  4,  solving  the  above  moment  problems,  -«»  <  k  <  ■».  Moreover,  the  result  of 
Ingham  implies  as  explained  in  [12],  [26],  that  for  each  k 

c*\2  <  /T  |uk(t)|2dt  4  C2N2 


where 


N2  =  2 1 k !  ( | k ]  +  1)lvk0(T)|2 


k  -  ±1,12,... 


+  j,  Vi,i|vk,t(T)|2  +  j,  Anti,*lvk7*(T>r 

*  j,  U|k|,t|wk,t(T)|2  +  j,  y|k|,llwk,t<T>|: 


2  2 

For  k  =  0  we  must  add  I Cqq (T) I  +  f wQ0 (T ) |  .  Since 


T  2s  «  T 

/  /  |u(8, t ) |  d8dt  “  £  /  lu.( t ) I  dt 


(6.46) 


0  0 


k— «•  0 
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we  see  that  the  above  moment  problems,  equivalent  to  the  control  problem,  can  be  solved 
with  (6.46)  finite,  provided  that 


which  is  the  same  as  saying  that  the  norm  of  the  final  state  in  H  should  be  finite.  Me 
have,  then,  the  exact  controllability  result  that  any  H  state  may  be  controlled  to  any 
other  H  state  during  a  time  interval  of  length  T  >  4  with  the  control  configuration  we 
have  described  here.  As  discussed  in  connection  with  the  wave  equation  in  [FFl,  [GG] ,  one 
cannot  be  sure  that  the  state  of  the  system  remains  in  H  for  all  t  6  I0,T).  However,  in 


the  present  case  of  the  Maxwell  equations  one  can  show  that  these  states  do  lie  in 

h  -  «yd<n>. 


7.  CONCLUDING  REMARKS 


Tha  approximate  controllability  reaults  of  Section  5  would  appear  to  be  extendable  to 
domaine  other  than  rectangular  ones  but  the  precise  method  of  extension  remains  to  be 
worked  out.  Me  will  indicate  some  aspects  of  this  problem  which  are  clear  from  our  current 
work. 

First  of  all,  the  result  of  Section  5  is  almost  trivially  extended  to  the  case  where 
control  is  exercised  only  on  a  subset  {(*,y)10  <a<y<b<w),  b  >  a,  of 
(<»,y)|0  <  y  <  »).  The  only  change  is  that  the  interval  1+5<t«T-1-6  appearing 
in  (5.14)  and  subsequently  must  be  sx>dified  to  d  +  6<t<T-d-4  where 

d  -  inf  {  sup  { [ (w  -  o2  ♦  (n  -  y>2]1/2}}  . 
a<y<b  0<(<« 

o<n<* 

If  ♦<*«y,t)  S  (s.y.t)  i  0  for  S  <  t  <  T  -  S,  a  <  y  <  b,  the  Holmgren  theorem  will 
still  apply  to  show  that  +(x,y,t)  =  0,  (x,y)  e  R,  d+S<t<T-d-5.  After  that  the 

remainder  of  the  proof  is  the  samei  the  same  eigenfunctions  and  frequencies  must  be  dealt 
with,  the  functions  sin  Jy,  cos  jy  are  still  independent  on  a  <  y  <  b  if  b>a  and 
the  conditions 

A 

Vg ( x , y )  -  0,  <x,y)  €  B  -  {<w,y)|a  <  y  «  b} 

iv0 

ST  (w'y)  “  °'  »  <  y  <  . 

still  show  Vg(x,y)  SO  in  R. 

The  first  limitation  of  the  method  which  we  have  used  in  Section  5  lies  in  its 
dependence  on  the  construction  of  +(x,y,t)  as  a  linear  combination  of  partial  derivatives 
of  v  and  w.  It  is  necessary  to  have  a  solution  of  the  wave  equation  to  which  Holmgren's 
theorem  may  be  applied.  This  part  of  the  proof  can  still  be  used  for  non-rectangular 
domains  as  long  as  a  portion  of  the  boundary  on  which  control  is  applied  is  a  straight  line 
segment.  Assuming  the  segment  parallel  to  the  y-axis,  one  can  construct  t  by  the  formula 
(5.8)  again  and  show  that  t  «n<J  both  vanish  on  the  straight  line  segment  in 

question,  allowing  subsequent  application  of  the  Holmgren  theorem  to  show  ♦(x,y,t)  =  0 


for  (X/y)  e  R  and  t  In  some  interval  d  ♦  i  <  t  <  t  •  i  •  I,  with  d  depending  on  the 
geometry  of  R.  But  then  we  are  faced  with  a  second  limitation. 

The  second  limitation  of  the  method  which  we  have  used  lies  in  its  reliance  on  the 
specific  form  of  the  eigenfunctions  and  frequencies  to  pass  from  g(x,y,t)  3  0  to  the 
conclusion  that  both  v(x,y,t)  and  w(x,y,t)  are  likewise  identically  aero.  It  needs  to 
be  emphasised  that  no  local  analysis  will  suffice  here.  In  the  one  dimensional  case  (see 
our  remarks  at  the  end  of  Section  3 )  if  the  control  problem  is  stated  for  boundary 
conditions 

v(0,t)  -  0,  ||  (1,t)  -  eni(t)  (7.1) 


§7  (0,t)  -  0,  |2  <1,t>  -  0u(t> 


(7.2) 


the  v,  w  constructed  as  in  Section  4  will  satisfy  the  wave  equation  and 


v(0,t>  -  0,  (1,t>  -  0  , 


(7.3) 


—■  U>,t)  •  0,  C1,t>  -  0  , 

ox  9x 


(7.4) 


a  (1,t)  +  0  jj'  (1»t)  =  ♦  ( 1 , t )  ■  0 


(7.5) 


Here  if  we  take  w  to  be  a  non-zero  solution  of  the  wave  equation  satisfying  (7.4)  and 
take 

v(x,t)  -  -  |  /X  <C.t)dC 

°  o  8t 

we  clearly  have  v(0,t)  “  0, 


57  (I,t)  -  |  j'  ~  (C,t)dC 
at  *  o  at2 


*  -  5  /’  *5  -  |  (1=  <o,t)  -  £  d.t>)  -  o  . 


ao  »52 


a  '8x 
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(x,t)  ■  - 


m  K  ,3“ 

|  /  M  <e,t)<* 

“  0  8t3 


a  X  .3* 

X  /  -^2  <«.t)dt 

a  o  itH 


so  that  v  aatlafiaa  the  wave  equation  and,  clearly,  (7.5)  is  also  satisfied.  Thus  the 
wave  equation  with  (7.1),  (7.2)  is  not  approximately  controllable!  +(x,t)  = 
a  (x,t)  ♦  B  (x,t)  S  0  but  this  does  not  imply  that  v  or  w  are  identically  equal 
to  zero.  The  additional  condition  which  makes  this  work  in  (3.7)  ff.  is  the  fact  that  one 


can  show  there  that 


-a  ( 0 , t )  ♦  B  (0,t)  •  0  . 


It  seems  likely  that  the  question  of  whether  or  not  4  ■  0  implies  that  both  v  and  w, 
equivalently  v  and  w,  are  both  zero  must  eventually  reduce  to  a  boundary  value  problem 
of  an  as  yet  unidentified  type. 

At  the  present  writing  there  is  only  one,  rather  curious,  result  which  we  can  offer 
which  yields  approximate  controllability  for  a  domain  R  of  rather  general  shape.  We 
suppose  that  the  "control  boundary"  Bj  B  “  3R  includes  two  nonparallel  line  segments, 
l 1  and  t2,  with  unit  exterior  normals  v,  and  Vj.  Proceeding  as  before  we  can  show, 
applying  the  Holmgren  theorem  together  with 


jit  "  °  °n  *1'  *2 


0,  i  »  1,2  on  respectively. 


a  ”  °«  i  *  1,2  on  respectively. 


that  both 


3v  -  3w 
♦i  +  b  it ' 


'*■  •**-**,-  .  ">+*<}■• - -  '»  ■  ■  ■ 


♦2 


*v_  .  3w 

*  e  H 


(7.7) 


must  vanish  identically  in  R  for  d+!<t<T'd-!,  4  >  0  arbitrary,  d  >  0 

depending  on  the  geometry  of  R  and  B,  the  location  of  1 1  and  t2  within  B,  etc. 

But  then  both  ^  and  $2  must  vanish  on  t 1  (say)  for  these  values  of  t.  Subtracting 
(7.6)  from  (7.7)  we  see  that 

*  A 

“(tir-  "  ■ 0  on  x  td  +  4,  t  -  d  -  4] 

12 

This  shows,  since  f1  and  i2  are  not  parallel,  that  a  nontangential  derivative  of  v 
vanishes  on  l.)x[d  +  6,T-d-4].  Combining  this  with  “  0  on  i1  and  applying 

the  Holmgren  theorem  to  v  alone,  much  as  in  [5],  (13),  we  are  able  to  conclude  v  =  0, 
provided  T  is  appropriately  large.  Then  one  easily  has  the  same  --  .-'it  for  w  and 
approximate  controllability  follows. 

This  result  gives  approximate  controllability  for  R  equal  to  the  interior  of  any 
closed  polyhedron  in  R  with  control  on  at  least  two  sides. 

Further  inspection  of  this  argument  shows  that  only  t2  needs  to  be  assumed  to  be  a 
line  segment.  That  is  needed  in  order  to  identify  42  as  a  solution  of  the  wave 
equation.  We  may  then  take  1 1  to  be  any  smooth  portion  of  which  is  never  parallel 

to  £ j  and  achieve  the  same  result. 

Finally,  let  us  indicate  that  we  are  very  much  aware  of  the  limitations,  from  the 
point  of  view  of  actual  implementation,  of  the  control  configuration  discussed  in  this 
paper.  In  principle,  at  least,  the  boundary  conditions  (1.7),  (1.8),  along  with  the 
further  "single  layer"  condition  discussed  in  connection  with  Figure  3.1,  could  be  achieved 
with  conducting  bars  attached  to  terminals  as  shown  in  Figure  3. 
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Figure  3.  Conducting  Bar  and  Busses 

The  perfectly  conducting  busses  perpendicular  to  the  boundary  of  f!  ensure  that  the  normal 
component  of  S .  Ev>  is  zero  just  outside  0,  provided  that  no  net  change  is  allowed  to 
accumulate  at  the  boundary  of  fl«  i.e.,  in  the  conducting  bar.  Thus  the  potentials  at 
C  and  D  must  be  regulated  so  that  the  potential  difference  C  -  D  ensures  the  correct 
controlling  current  through  the  surface  bar  B  while  C  +  D  is  set  so  that  there  is  no 
accumulation  of  charge  at  the  bounding  surface. 

We  have  not  considered  any  effects  of  propagation  delays  in  the  controlling  circuits  - 
i.e.,  we  have  not  assumed  that  these  are  distributed  parameter  systems.  This  assumption, 
and  evident  limitations  on  the  speed  with  which  prescribed  currents  can  be  computed  and 
established  in  the  controlling  circuits  together  with  sensing  limitations,  place  admittedly 
severe  limitations  on  what  can  be  done  "open  loop".  It  Is  likely  that  the  eventual 
significance  of  our  results  will  be  most  evident  in  connection  with  closed  loop  behavior 
wherein  time  varying  mar  letic  fields  8  near  the  boundary  of  8  induce  currents  in  the 
bars  B  which,  being  resistive,  will  then  act  as  energy  dissipators.  we  hope  to  discuss 
this  topic  in  later  work. 
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Another  control  configuration  nay  be  obtained  by  supposing  the  boundary  of  A  to  be  a 


perfectly  conducting  sheet  of  material  to  which  electromagnets  are  attached  in  a  dense 
array  as  shown  in  Figure  4. 


Figure  4.  Electromagnet  Array 


If  J  denotes  the  current  through  the  windings  of  the  electromagnets,  then  we  shall  have 


0 


and 


Hv  -  oJ 

where  a  is  dependent  on  the  electromagnet's  configuration.  The  theory  in  this  case  will 


take  much  the  same  form  as  the  one  discussed  in  this  paper 
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